Observation of the ground-state-geometric phase in a Heisenberg XY model 
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Geometric phases play a central role in a variety of quantum phenomena, especially in condensed 
matter physics. Recently, it was shown that this fundamental concept exhibits a connection to 
quantum phase transitions where the system undergoes a qualitative change in the ground state 
when a control parameter in its Hamiltonian is varied. Here we report the first experimental study 
using the geometric phase as a topological test of quantum transitions of the ground state in a 
Heisenberg XY spin model. Using NMR interferometry, we measure the geometric phase for different 
adiabatic circuits that do not pass through points of degeneracy. 
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When a quantum system is subjected to a cyclic adi- 
abatic evolution, it returns to its original state but may 
acquire a geometric phase factor in addition to the dy- 
namical one. Berry made this surprising discovery in 
1984 [1 , so that this is also known as Berry's phase. 
Later this phase was generalized in various directions, to 
include a more general case of noncyclic and nonadiabatic 
evolution [2], and even the case of mixed states. Geomet- 
ric phases (GP) have been observed in a wide variety of 
physical systems, e.g., in spin-polarized neutrons [3], nu- 
clear magnetic resonance (NMR) [4 and superconduct- 
ing systems [3]. Moreover, GP has found applications 
to many areas, such as molecular dynamics, many-body 
systems and quantum computation [B] [7] . 

Very recently, the GP of many-body systems has been 
shown to be closely connected to quantum phase transi- 
tions (QPTs), an important phenomenon in condensed 
matter physics [8] [9]. QPTs occur at zero tempera- 
ture and describe abrupt changes in the properties of the 
ground state resulting from the presence of level crossings 
or avoided crossings [10] . Recently, different methods re- 
lated to quantum information have been developed for 
characterizing QPTs, including the fidelity [11], quantum 
entanglement [I2j [13] and some other geometric proper- 
ties [14]. The GP, which is a measure of the curvature 
of Hilbert space, can reflect the energy level structure 
to fingerprint certain features of QPTs. Carollo and Pa- 
chos [8] demonstrated that the GP difference between 
the ground state and the first excited state encounters 
a singularity when the system undergoes a QPT in the 
XY spin chain. Zhu [9] revealed that GP associated with 
its ground state exhibits universality, or scaling behavior, 
around the critical point. Besides the study in the ther- 
modynamical limit, it was also shown that the GP could 
be used to detect level crossings for a two-qubit system 
with XY interaction [15] . As a complement to these theo- 
retical investigations, it appears highly desirable to have 
experimental evidence for these effects. 

In this Letter, we report the first experiment that 
shows this important connection between the GP and 



the energy level structure (i.e., level crossing points) in 
a Heisenberg XY spin model. In our experiment, the 
system Hamiltonian changes adiabatically along a closed 
trajectory in parameter space while the system, which 
is in the ground state of the Hamiltonian, accumulates 
a GP. Depending on the region in parameter space, the 
resulting GP is zero or has a finite value. These regions 
in parameter space are separated by a line where the 
ground state of the system becomes degenerate [15] . Us- 
ing adiabatic state preparation and NMR interferometry, 
we observe the transitions of GP on both sides of the 
level crossings point. This experiment might be viewed 
as a first meaningful step to use GP as a fingerprint for 
ovserving QPTs. 

Consider a one-dimensional spin- 1/2 XY model in a 
uniform external magnetic field along the z axis: 

= -£(^m +1 + ^^ +1 )-^£< 
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where (J^iy = x, y, z) denote the Pauli matrices for qubit 
A is the strength of the external magnetic field, and 
7 measures the anisotropy of the coupling strength in 
the XY plane. This model is exactly solvable and can 
be diagonalized by the Jordan- Wigner transformation, 
Fourier transformation and then Bogoliubov transforma- 
tion [IB] . However, it still contains a rich phase structure 
[TO] , Barouch and McCoy [17] investigated the statisti- 
cal mechanics of this model in the thermodynamical limit 
and showed that a circle (A 2 + 7 2 = 1) separates the os- 
cillatory phase (inside) from the para- or ferro-magnetic 
phase (outside) . At the level crossing or avoided crossing 
between ground state and first excited state, the ground 
state changes discontinuously. As a result, the GP associ- 
ated with the ground state also changes discontinuously. 
Theoretical work has demonstrated the close relation be- 
tween GPs and the energy level structures, thereby re- 
vealing the ground-state properties [8] [9], even in the 
two-qubit case [IB] . 

We now consider the GP that results in this system if 
the Hamiltonian rotates around the z-axis, H(A,7,</>) = 
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U$(<I>)H(\, -y)U z (</>) with [^(0) = Uk e ~^ akz [8 • # has 
the same spectrum as independent of 0. Here we 
study a minimal model of two qubits coupled by an XY- 
type interaction [15]. The eigenvalues of ti are ±1 and 
±r, where r = y^A 2 + 7 2 . The ground state is 
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(|01) + |10)), r<l 



(1) 



cos-|00) +sin-e- i2 ^|ll), r > 1 



where tan# = 7/ A. For r < 1, the ground state is thus 
invariant; for r = 1, it is doubly degenerate; and for 
r > 1, it is spanned by the two states 1 00) and |11), with 
coefficients that depend on the angle 6. 

If we let the Hamiltonian travel along a cyclic 
path in the parameter space (A, 7,0), we can consider 
the subspace spanned by 1 00) and |11), which con- 
tains the ground state, as a pseudo-spin 1/2, where 
the spin evolves in an effective magnetic field B = 
r (sin cos 20, sin 6 sin 20, cos 6) . Using the standard for- 
mula f3 g = i (\£ g \d(f)\^ g ) [1 , the ground state accumu- 
lates a GP 



^(0:O^tt) 



0, r<l, 
7r(l — cos #), r > 1. 
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As shown in Fig.[TJa), it is useful to represent the trajec- 
tory in a parameter space spanned by 7 cos(20), 7 sin(20) 
and A. Here, the sphere with radius r = 1 marks the 
points where the Hamiltonian is degenerate. Inside this 
sphere (r < 1), the GP vanishes, while it has a finite 
value that depends on the opening angle of the cone 
subtended by the circuit. A special case is the XX spin 
model (i.e., 7 = 0). Here, the GP always vanishes, be- 
cause the operation U z does not change the Hamiltonian 
of the system. While we are considering here only a min- 
imal two-spin model, the ground state and the ground 
state energy of the XY model in the thermodynamic limit 
are similar [T7] . 

When the system undergoes the cyclic adiabatic evo- 
lution along there will also be an additional dynamic 
phase generated, relative to the instantaneous energy of 
the system, besides the GP. Hence, in order to acquire 
the pure geometric part, we have to eliminate the dy- 
namical contribution. To eliminate the dynamical con- 
tribution to the phase shift, we combine two experi- 
ments with the closed paths C and C [2], which gener- 
ate the same geometrical phases, but opposite dynam- 
ical phases. The two trajectories have the same geo- 
metrical shape (cones), but their Hamiltonians T~L and 
— H and thus their dynamical phases add to zero. Dur- 
ing the first period, the Hamiltonian H(A,7, 0) follows 
the closed curve C in the parameter space r = (r, 0,0), 
with changing from to 7r, as schematically shown 
by the the red circles (labeled by C in the upper part) 
for A > in FIG. [lja). During the second period, the 
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FIG. 1: (Color online) (a) Parameter-space representation of 
the cyclic adiabatic evolutions that generate the GP. Two 
closed paths C and C, related by inversion symmetry, were 
combined for observing a purely GP. The cycles are horizonal, 
i.e., A is constant and 7 is constant. The observed GP depends 
on the angle if the circles are outside the sphere r — 1 (shown 
in black) and vanish if the curves are inside the sphere, (b) 
Energy level diagram of the time-dependent H a d(s) for ASP 
(denoted by the dashed lines), and the optimal function of adi- 
abatic parameter s(t) (denoted by the solid line) calculated 
for a constant adiabaticity factor, when 7 = 0.5, A = 0.9920. 
The black dots represent the experimental values for the dis- 
cretized scan. 



Hamiltonian — ti = R^ kz (7r)H(—X, 7, (j))Rk z {^) follows the 
curve (7, shown in the lower part of FIG. Pa a). Here 
Rkz(n) — e~ q, 2 a z^k = 1 or 2) rotates one of the two 
spins around the z-axis. For the circuit C, the result- 
ing phase is j3c = 7r(l — cos#) — rT, where T is the 
cycle time, where we have assumed r > 1. For (7, 
Pc ~ ^(l — cos#) + rT because the sign of the eigen- 
value of the state \^ g ) changes for —T~L. The sum of the 
two phases, fic + Pc = 2^r(l — cos#) is thus purely geo- 
metrical. If r < 1, the dynamical component changes to 
— T for C and T for C while the GP vanishes. 

For measuring the GP, we use NMR interferometry 
[4j [18]. This requires an ancilla qubit that is coupled 
to the system undergoing the circuit. FIG. [2] shows 
schematically the experiment, including the adiabatic 
state preparation (ASP) of the two qubit system into the 
ground state of the Heisenberg XY model, and the gen- 
eration of a superposition of the ancilla qubit by a Had- 
mard gate. The subsequent adiabatic circuit L^, which is 
conditional on the state of the ancilla qubit, implements 
the interferometer^ = |0)(0| a (8)l + |l)(l| a (8)^, where 1 
represents a 4 x 4 unit operator and the unitary operator 
Ui is the cyclic adiabatic evolution on the system qubits 
along the chosen path C or C. The phase acquired dur- 
ing this path appears then directly as a relative phase 
in the superposition of the two ancilla states and can be 
measured in the NMR spectrum of the ancilla spin. 

The experiment was carried out on a Bruker Avance III 
400MHz (9.4 T) spectrometer at the room temperature. 
The three qubits 0, 1 and 2 in the quantum circuit (FIG. 
[2^a)) were represented by the X H, 13 C, and 19 F nuclear 
spins in Diethyl-fluoromalonate. The relaxation times for 
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FIG. 2: Quantum circuit for measuring the ground- state GP. 
H is the Hadamard gate, and following the adiabatic state 
preparation (ASP), the operation Ui performs a cyclic adia- 
batic evolution of the system qubits (1 and 2) conditionally 
when the ancilla path qubit is in the state |1). 



Is. The natural Hamiltonian of 
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all three spins are T2 
this system is Unmr = ~ J2i=o ^f a \ 
where uji is the Larmor frequency for spin i and J 2 
the coupling constants J01 = 160.7 Hz, J12 
Hz and J02 = 47.6 Hz. As the sample is not labeled, 
the relative phase information on 1 H at the end of the 
quantum circuit was obtained through the 13 C spectrum 
by a SWAP operation between 13 C and 1 H [T3] . 

In the experiment, we first initialized the system into 
the pseudopure state (PPS) p oo = ^pl + e|000)(000| 
by spatial averaging [13], with the polarization e « 10 -5 . 
Then we prepared the ground state of the Heisenberg XY 
Hamiltonian by an adiabatic passage: A rf pulse rotated 
the spins from the z- to the — x-axis, i.e., to the ground 
state of Ho = a x> and then this Hamiltonian was 
slowly changed into the target XY Hamiltonian H (A, 7), 
always fulfilling the adiabatic condition k <C 1 [19 . This 
assures that the resulting final state is close to the desired 
ground state of the XY model. We optimized the time 
dependence of the transfer by choosing H a d(t) = [1 — 
s(t)]H + s(t)H with < s(t) < 1. The solid line in 
FIG.[TJb) shows the corresponding time dependence for a 
constant n. The time dependence of s(t) was chosen such 
that the adiabaticity parameter n < 0.25 at all times. 

In the experiment, the adiabatic transfer was per- 
formed in discrete steps. The parameter s(t) there- 
fore assumes discrete values s m with m = 0, ...,Mp, 
and for each period of duration (5, the correspond- 
ing Hamiltonian 7-L a d[sm] was generated by a mul- 
tiple pulse sequence: Up (5) = e -^ad[s m )s _ 



0(£ 3 ), 
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use of Trotter's formula. S and Mp were chosen by simul- 
taneously considering this stepwise approximation and 
the adiabaticity criterion. The experimental values s m 
for the discretized scan are represented by black dots in 
FIG. 1(b). The theoretical fidelity of this stepwise trans- 
fer process was > 0.99, and the experimental fidelity was 
> 0.98. 

After the preparation of the ground state, we applied 
the cyclic adiabatic variation C or C. The corresponding 
control operation tic or Uq was generated in the form 
of a discretized adiabatic scan, as described for the ASP 
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FIG. 3: (a) Experimental NMR spectra for two specific pa- 
rameter sets, Pi for 7/(0,0.5) and P 2 for H(0.878, 0.5). From 
top to bottom, the spectrum corresponds to the initial ground 
state |^ P (A,7)), Uc\^g(\j)) for the adiabatic path C and 
Uc\^ g (\j)) for the adiabatic path C. (b) Measured ground- 
state GP of the Heisenberg XY model (points) for different 
parameter sets (A, 7) compared to the theoretical expectations 
(solid curves). 



part. Again, the parameters of the scan were optimized 
to keep the fidelity> 0.99. At the end of the scan, the 
accumulated phase was measured. 

Fig. [3] (a) shows two representative examples of the 
resulting data: The spectra on the left hand side cor- 
respond to the states before the adiabatic circuit, after 
traversing the circuit C, and after traversing C for the 
Hamiltonian parameters (A, 7) = (0,0.5). Clearly, in this 
situation, we are within the sphere r = 1, so the ground 
state of the system is A= (| 01) + 1 10) ) . This is verified 
by the experimental data, where only the two resonance 
lines are visible that correspond to the states 1 01) and 1 10) 
of the system. In the initial state, the two lines appear 
in absorption; this corresponds to the reference phase 
(p = 0. During the circuit C or (7, which is traversed 
over a time T = 3, the system should acquire a phase 
T. In the experimental data, we find that the lines are 
inverted; a numerical analysis of the lines yields phases of 
Pt{C) w 170.6° and -173.0°. Thus the result- 

ing GP p g = (/3 t (C) + Pt{C))/2 « -1.2°, which is close 
to the theoretically expected value of zero. The right- 
hand part of Fig. [3ja) shows the corresponding results 
for an adiabatic circuit outside of the sphere r = 1, where 
we expect to observe a non-vanishing GP. In this case, 
we observe clearly different phases for the two circuits, 
whose duration is now T = 10.7. The measured phases 
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are f3 t (C) « -62.6° and (3 t (C) « 114.2°, corresponding 
to a GP of 25.8°. 

Fig. [3] (b) shows the GP measured for different param- 
eters (A, 7). The symbols show experimental data points, 
while the curves that connect the points show the theo- 
retically expected GP as a function of the magnetic field 
strength A, for a constant anisotropy parameter 7. In all 
cases, the observed GPs are compatible with the theoreti- 
cally expected values: zero if the parameters (A, 7, cp) fall 
inside the sphere with radius r = 1, a sudden increase 
to the maximum value just outside the sphere, where 
the opening angle of the cone subtended by the circuit 
reaches a maximum, and then decreasing as the circuit 
C is moved away from the origin. Increasing values of 7 
correspond to larger circles C and thus bigger values of 0. 
The points marked P\ and P2 correspond to the spectra 
shown in the upper part of the figure. 

The relevant sources of experimental errors mainly 
came from undesired transitions induced by the time- 
dependent Hamiltonian, inhomogeneities of rf fields and 
static magnetic fields, imperfect calibration of rotations 
and relaxation. We used a numerical optimization pro- 
cedure to minimize undesired transitions during the adi- 
abatic passage. The durations of individual experiments 
ranged from 30 ms to 90 ms, short compared to the re- 
laxation time T2 ~ Is. The experimental error of the 
geometric phase was less than 3°. The imperfection of 
the initial state would also contribute to this. Using the 
experimentally reconstructed density matrices for the ini- 
tial states, we found that this effect contributed « 1° to 
the errors. 

In summary, we have detected the ground-state GP 
in the Heisenberg XY model, after preparing the ini- 
tial state by an adiabatic passage. The Heisenberg XY 
model was simulated by a multiple-pulse sequence, and 
the phase was measured by NMR interferometry. Our 
proof-of-principle experiment illustrates that the ground- 
state GP can serve as a fingerprint of the energy-level 
crossing points that result in a QPT in the thermody- 
namic limit. The ground- state GP is a robust indicator 
that is immune to some experimental imperfections [20] 
and provides an experimental method that does not need 
to cross the critical point. 

It would be very interesting to extend this experiment 
to larger spin systems. For this, two issues are relevant: 
(i) the effectiveness of the ASP and {%%) the realization of 
quantum circuit consisting of a quantum interferometer 
and quantum simulation. For the first issue, although a 
decisive mathematical analysis of the efficiency of ASP 
is difficult, numerical simulations (up to 128 qubits) [21] 
indicate a polynomial growth of the median runtime of 
an adiabatic evolution with the system size. On the sec- 
ond issue, quantum interferometry has become a mature 
technique, and the Heisenberg XY model has been ef- 
ficiently simulated by a universal quantum circuit only 
involving the realizable single- and two-qubit logic gates 



[22] , Moreover, the diagonalization theory of the XY 
model shows a valid energy gap between the two lowest 
energies, which guarantees the viability of the cyclic adi- 
abatic evolution to generate the ground-state GP, even in 
the thermodynamic limit [8] . Recent research also shows 
that a 10-qubit system already represents a good approx- 
imation to the thermodynamical limit [23]. Therefore, 
the present scheme is in principle applicable to larger 
spin systems, when the technical difficulties in building 
a medium-scale quantum computer are overcome. This 
significant connection between GPs and QPTs is not a 
specific feature of the XY model, but remains valid in 
a general case [8j [9]. We hope that this experimental 
work will contribute to an improved understanding of the 
ground-state properties and QPTs in many-body quan- 
tum systems. 
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FIG. 4: (Color online) Relevant properties of the quantum 
register, (a) The molecular structure of Diethylfluoroma- 
lonate and the oval marks the three spins used as qubits. 
(b) The relevant NMR parameters: the resonance frequencies 
uji (on the diagonal), the J-coupling constants Jij (below the 
diagonal), and the relaxation times T\ and T<i in the last two 
columns, (c) NMR spectrum of the 13 C obtained through a 
read-out tt/2 pulse on the equilibrium state, where the four 
resonance lines are labeled by the corresponding states of the 
two other qubits. 
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FIG. 5: Pulse sequence for implementing the SWAP opera- 
tion on spins en and aj with the corresponding matrix rep- 
resentative. Narrow rectangles indicate tt/2 pulses and wide 
rectangles tt pulses, and the pulse phases are indicated above 
them. 



Experimental procedure 

The experiment, summarized in Fig. 2(a) in the pa- 
per, includes three steps: (i) Adiabatic state preparation 
(ASP): to prepare the ground state of the XY spin model 
H(A, 7) at different points of the Hilbert parameter space 
(A, 7) by adiabatic evolution; (ii) NMR interferometry: 
to generate the pure geometric phase (GP) on one of the 
two paths where an auxiliary spin is introduced; (Hi) 
Phase measurement: to obtain the GP by measuring the 
relative phase between the two paths via quadrature de- 
tection in NMR. 



1. Adiabatic state preparation 



SUPPORTING ONLINE MATERIAL 
Quantum simulator and characterization 

For the quantum register for these experiments, we 
selected the 1 H, 13 C, and 19 F nuclear spins of Diethyl- 
fluoromalonate dissolved in d-chloroform. The relevant 
system parameters are listed in Fig. [4] Experiments were 
carried out at room temperature, using a Bruker Avance 
III 400 MHz (9.4 T) spectrometer equipped with a QXI 
probe with pulsed field gradient. 

Because we used an unlabeled sample, the molecules 
with a 13 C nucleus, which we used as the quantum regis- 
ter, were present at a concentration of about 1%. The 1 H 
and 19 F spectra were dominated by signals from the 2- 
qubit molecules containing the 12 C isotope, while the sig- 
nals from the quantum register with the 13 C nucleus ap- 
peared only as small (0.5%) satellites. To effectively sep- 
arate this signal from that of the dominant background, 
transfered the state of the 1 H and 19 F qubits to the 13 C 
qubit by a SWAP gate and read the state through the 
13 C spectrum. The matrix representation of the SWAP 
operation for two spins ai and aj is shown in Fig. |5j 
together with the corresponding pulse sequence. 



To prepare the system in the ground state of a Hamilto- 
nian H(X, 7), the system can start with an initial Hamil- 
tonian H(0) = Ho, whose ground state \ip g (0)) is known 
and then it is adiabatically driven to the target Hamil- 
tonian H(T) = %(A, 7). A time-dependent Hamiltonian 
H a d(t) smoothly interpolates between H(0) and H(T): 



H ad (t) = [1 - s(t)]H + s(t)H(\,j) , 



(3) 



where the function s(t) varies from to 1 to parametrize 
the interpolation. If the quantum system starts in |^(0)) 
and the variation of T~L(t) is adiabatic, the final state 
reached will be close to the ground state of H(T) = 
H(A,7). To ensure that the system is prepared in the 
ground state of %(A, 7), the sufficiently slow variation of 
H a d(t) means that the traditional adiabatic condition [? 



e le (t)-s g (t) 



< 1 



(4) 



is fulfilled, where \^ g (t)) and \^i e {t)) refer to the instan- 
taneous ground state and the first excited state, respec- 
tively, and £ g (t), Si e (t) are the corresponding energies. 

To do this, we chose the initial Hamiltonian of the 
system = <j 2 + <r 3 , whose ground state 1^(0)) 23 = 



71 (I °) - I 1 )) ® 71 (1°) - I 1 )) is well-known. Starting from 
the thermal equilibrium state, the system was first ini- 
tialized into a pseudopure state (PPS) pooo = + 
e 1 000) (000 1 by spatial averaging [25], where 1 represents 
the unity operator and e « 10 -5 represents the thermal 
polarization. The |0) and |1) states correspond to the 
two eigenstates of o z the spin- up and spin-down states, 
respectively. The normalized deviation density matrix of 
the PPS [25 , pa = [p— (1 — was reconstructed by 

quantum state tomography [25j [26] , which involves the 
application of 7 readout pulses and recording of the spec- 
tra of all three channels to obtain the coefficients for the 
64 operators comprising a complete operator basis of the 
three-spin system. The experimentally determined state 

fidelity was F = , T ^£z* Pth \ 2 , « 0.99. Then the 

initial state |0)i (g> \ip g (0)) 23 was prepared by two pseudo- 
inverse-Hamdamard gates, i.e., [7r/2]- y pulses on the sys- 
tem qubits 2 and 3. 

An optimal function of s(t) determines the efficiency of 
ASP. To find the optimal interpolation function s(t) for 
the adiabatic process, we rewrite the adiabatic condition 
of Eq. Q as 



ds(t) 



dt 



\eie{t) - e g {t)\* 



dU ad {s) 
ds 



IV>le(*)> 



X, 



(5) 



which defines the optimal sweep of the control param- 
eter s(t) with the scan speed The required time 
dependence of s(t) was numerically optimized for con- 
stant adiabaticity parameter k = represented by 
the solid line in Fig. 1 (b) in the paper for the Hamil- 
tonian H(\,j) = H(0.992, 0.5). The time dependence 
of s(t) was chosen such that the adiabaticity parameter 
k < 0.25 at all times. 

For the experimental implementation, we have to dis- 
cretize the continuous adiabatic passage into Mp + 1 seg- 
ments (i.e., s m = s(-^Tp) with m = 0, ...,Mp, so = 
and sm p = 1), and generate the instantaneous discretized 
Hamiltonian H a d[sm] = [1 — s m ]Ho+s m H(A, 7) for a time 
S. The evolution operator for the rath step is given by 



U„ 



i5U a 



where 5 = T P /(M P + 1). The total evolution is 



M P 



Up = H U n 



(6) 



(?) 



ra=0 



Since Hq and H(\,j) in 7-L a d[sm] do not commute, the 
operator U m (S) is approximately implemented by the use 
of Trotter's formula [24] : 

U m e -iSH ad [Sm] = e -i^[lSm]Ho 

xe -iSs m n(x,j) e -i^ [i-s m ]n _^ o(5 3 ). (8) 



h_ ■ ■ m ■-- y - x — 

f_r i__J___i 1. 



71 pulse 



^7l/2 pulse [u] r| puis 



FIG. 6: The pulse sequence for a step of adiabatic state prepa- 
ration. Here, r± l^rl, V = (1 ~ s m )S, r}- = f - 77, 
a - (-r - l)s m £, /? - (-r + l)s m 5. 



For this stepwise approximation, the duration of each 
time step S has to be chosen such that (i) the time S 
is short enough that the Hamiltonian simulation holds 
and (ii) the adiabaticity criterion remains valid, i.e., the 
total time Tp is long enough. The experimental values 
Sm for the discretized scan are represented by black dots 
in FIG. 2(b) in the paper, which keeps a high theoret- 
ical fidelity (more than 0.99) of the final state for this 
stepwise transfer process of our ASP. The two opera- 
tions e - *2 [ 1-s ™]^o anc [ e -i5s m H(\,f) can fog precisely sim- 
ulated: e -z 2[ 1-s ™]^° can be easily realized using NMR 
radiofrequency pulses, while we use quantum techniques 
to simulate the XY Hamiltonian: 



Vie 



■i8smHd(X,^) 



v d , 



where the diagonal matrix %d{\l) = VdT-i(^,7)Vj: 



( -\A 2 +7 2 \ 
-10 
1 

V 00 v/P + f J 



^AHf + 1 ! x/^ + ^-l 2 



^(9) 



with 



v d 



( cos((9/2) - sin(0/2) \ 

1/V2 -1/V2 

1/V2 1/V2 

V sin(0/2) cos((9/2) / 



■ 6 + 7T/2 1 2 • 0-tt/2 1 2 



(10) 



Here tan# = 7/ A. Thus the operator U m can be imple- 
mented using a multi-pulse sequence [27], shown in Fig. [6] 
. In the case of the unsolved model for the target Hamil- 
tonian, the propagator still can be obtained by average 
Hamiltonian theory [28] . 

To confirm the success of ASP, we performed quantum 
state tomography on the final state at the end of adia- 
batic passage. For example, the experimentally recon- 
structed deviation density matrices of the system at the 
positions P 1 (7 = 0.5, A = 0) and P 2 (7 = 0.5, A = 0.878) 



7 




Real part 



Imaginary part 



FIG. 7: Experimental and theoretical deviation density ma- 
trices of the system after ASP at the positions Pi (7 = 0.5, A = 
0.878) and ^2(7 = 0.5, A = 0, along with theoretical predic- 
tions). The left and right columns denote the real and imag- 
inary components, respectively. 



are shown in Fig. [7| along with the theoretical expecta- 
tions. From the tomographically reconstructed density 
operators, we determined the experimental fidelity for 
our prepared states: = 0.98, Fp* p = 0.99. This 

proves that ASP prepared successfully the ground state 
of the XY model. 



2. NMR interferometry: A purely GP generation by 
conditionally cyclic adiabatic variations of the 
Hamiltonians 



Cyclic adiabatic evolution 



by 



The studied physical system is described 

n\i,f) k = ut{4>)u(\i)u z {4>) with u z (4>) = 

Y\ k e ~ l 2 a *, which has the same spectrum as %(A, 7), 
independent of 0. Its matrix form is 



H(A, 7 , 



A 7e i2< ^ 
10 
10 
-A 



\ 



(11) 



Clearly we can consider the subspace spanned by 1 00) 
and 1 11), which contains the ground state, as a pseudo- 
spin 1/2, where the spin evolves in an effective mag- 
netic field B = r(sin#cos20, sin # sin 20, cos#) with 
tan0 = 7/ A, while the subspace spanned by 1 01) and 
1 10) is independent of 0. Here, we adiabatically vary 
the angle 0, that is, the Hamiltonian depends on time 
through a set of parameters R(t) = R(#,0(t)) = 
(sin#cos20(t), sin#sin20(£), cos#). Thus we are inter- 
ested in the adiabatic evolution of the system as R(t) 



moves slowly along a path C in the parameter space, i.e., 

C : [0,T] S 2 with points R(t) e C, 

where the unit sphere S 2 = {R E 7Z 3 : |R| = 1} is 
the parameter space of the system. The quantum adia- 
batic theorem predicts that a system initially in one of its 
eigenstates |n; R(0)) will remain its instantaneous eigen- 
state |n;R(t)) of the Hamiltonian H(R(t)) throughout 
the process. 

If the parameters R(t) adiabatically traverse a closed 
path C : [0(0) = 0, 4>{T) = tt] and return, after some 
period T, to their original values: 



C : R(0) R(t) -> R(T) = R(0), 



then 



H(H(T)) = W(R(0)) 
E n (R(T)) = E n (R(p)) 
|n;R(T))(n;R(T)| = |n; R(0))(n; R(0)|. 

Here the Hamiltonian H(R) has the spectral resolution 
H(TL(t)) = ^£7 n (R(t))|n;R(t))(n;R(t)|. 

n 

However, the basis vectors |n;R) themselves in general 
not be unique over the whole parameter space. A new set 
of eigenvectors |n;R)' can be obtained by gauge trans- 
formations: 

|n;R)' = e^ R V;R), 

where Cn(R) are arbitrary real phase angles. One can use 
different parameterizations over different patches of the 
parameter space. Here we require that the closed path 
C is placed into one single patch O C M and the basis 
functions |n;R) is smooth and single- valued. 

Consequently, the cyclic adiabatic evolution along the 
closed path C : [0(0) = 0, 0(T) = tt] let the Hamiltonian 
T-L(H(t)) and the adiabatically evolving state \i(j(t))(ijj(t)\ 
return to their original forms in the parameter space as 
time progresses from t = to the period t = T. If our 
system start with the ground state \^ g (0)) = cos 1 1 00) + 
sin 1 1 11) of 7tC(R(0)), the adiabatically evolving state at 
time t is 

= eW°+M\9 g ((l>)) = eW° + Mul((l>)\9g(0)). 

The cyclic adiabatic evolution along the closed path C 
will generate the dynamical phase 

P d =<f> E g (t)dt=J> -y / A 2 +7 2 ^ = -^V / A 2 +7 2 
Jo Jo 

and the geometric phase (or Berry phase) 

P g = i£(* g (<P)\d rp \* g (<P)) = 7T(1 - COsO). 
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Berry phase for the GP obtained in the adiabatic ap- 
proximation is associated with a closed curve in the 
Hamiltonian parameter space [29 . After the cyclic adi- 
abatic evolution along the closed path C, the system 
state obtains the total phase /3t(C) = fid(C) + /3 g (C) — 
7r(l — cos#) — TyJ A 2 + 7 2 , mixing the GP with the dy- 
namical phase. 



A purely GP generation: eliminating the dynamical phase 

In order to obtain a purely GP, we have to design a 
reverse dynamical process to cancel out the dynamical 
phase but double the geometric one [U? ]. This can be 
implemented by another closed path C ((/) : tt — > 2tt) 
along the Hamiltonian — H(A,7, </>). In this case, the ini- 
tial state \ty g ) = cos ||00) + sin §|11) is not the ground 
state of — 7t£(A,7, 0), but the eigenstate of the Hamilto- 
nian — ?i(A,7, <fi) with the highest eigenvalue + 7 2 . 
Therefore the adiabatic passage was performed on this 
eigenstate along — H(A,7, qb). Thus we have the dynam- 
ical phase 



Pd(C) = j^ v / A 2 +7 2 ^ = ^\/A 2 +7 2 , 
and the geometric phase 



P 9 (C) = ij> (%(<l>m\%(<l>)) = tt(1 - cos0), 
which leads to the total phase f3 t (C) — /3d(C) + /3 g (C) = 

tt(1 - COS 0) + Ty^X 2 + 7 2 . 

The resulting effect of the two closed paths C and C 

is 

Pt{C) + = 2^ = 2tt(1 - costf). 

Consequently, the dynamic phase vanishes and we obtain 
a purely GP. 



NMR interferometry 

The phases generated can be detected by NMR in- 
terferometry, which consists of a Hadamard gate and a 
controlled- Z7i operation (see Fig. 2 (a) in the paper). The 
Hadamard gate is represented by the Hadamard matrix: 



H 



i 

V2 \1 "I 



which maps the basis state |0) to 4^(|0) + |1)) and |1) 
to 775(10) — The auxiliary qubit a was put into a 

superposition state A= (|0) + |1)) from the |0) state by a 
pseudo-Hadamard gate [7r/2]- y . 



Then the system state adiabatically traces out a closed 
path C (qb : —> tt) along the Hamiltonian T~L{\ 7, 0), 
but only if the auxiliary qubit is in state when the 
auxiliary qubit is in state |0), the system state is not 
affected. This can be realized by a controlled- Uc opera- 
tion: tic = |0)(0| a (g) 1 + |l)(l| a (g) /7c, where 1 represents 
a 4 x 4 unit operator and the unitary operator Uc is the 
cyclic adiabatic evolution on the system qubits along the 
chosen path. It effectively introduces a relative phase 
shift between the initially prepared superposition with 
known phase of the states of the auxiliary qubit when the 
cyclic adiabatic evolution Uc creates a non-zero phase. 
The process of the interferometer can be described as 

H - S(|0) o + |l) o )|^) 12 



|0>a|* 9 >12 



Un 



V2 
1 

71 



ma + e ilMC) + MC)] ll)a)l y g) 



9/12- 



Likewise, for the closed path C, The process is 

|0> o |*p>12 ^ -4d°>a + |l>a)l* fl >12 



V2 
1 

71 



The resulting effect of these two experiments results in 
1 



|0)a|* fl >12 



(|0}a 



l>a)|* fl )l2- 



where the purely GP is obtained by summing the relative 
phases in these two experiments. 



Experimental implementation for conditionally cyclic 
adiabatic evolutions along C and C 

The Hamiltonian H(A,7,0) varies adiabatically along 
the trajectory C, i.e., <\> changes slowly from to tt. Like 
in APS, the continuous Hamiltonian H(A,7, <j)) is dis- 
cretized into M + 1 steps in the range of <\> : — >• tt 
in the actual implementation. Likewise, we numerically 
optimized the adiabatic steps M + 1 and evolution time 
T to achieve a high fidelity of the instantaneous state of 
the system. In experiment, we chose M < 6 and T ~ 10 
which results in a theoretical fidelity of > 0.99 for both 
trajectories C and C. 

The unitary operation of the rath adiabatic step for a 
constant (j) m can be realized by the following decomposi- 
tion: 

; 1 /1 o a 



Jjm _ e -ii(l a -<rJ)«.W(A, 7 ,0 m )T 



= i a ® ul(4> m )e- i ^ ia -^ H ^ T i a ® u z (4> m ) 

= l a ® t/t (0m) yj e -^(l«-^)^(A )7 )r 1 a VdUzi<f)m ), 

where r = T j (M + 1). The total evolution is 



M 



U C (T) =]JU ( 



c ' 



(12) 



ra=0 
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FIG. 8: Pulse sequence for implementing the control op- 
eration Uc for the adiabatic path C. Here r± = \ |, 

4 = 1^1 (* = l,2),0 ± = f ±<f>, and a± = ^ 2 . 



Ul((/)m) is a rotation of the system qubits 
around the z axis, which can be realized by 
UkR kx (7r/2)R ky ^ m )R kx (-7r/2) and V ] d by Eq. 
(TlOj) . The conditional operation 

|(l°-^)®Kd(A,7)r 



= e 



e 62 



is also implemented by rf pulses and J-coupling evolu- 
tions, shown in Fig. [8] 

For the trajectory C, the propagator is 



c 



-i^(l a -a^)(g)[-'H(A,7,0^)]r 
-*5(1 O -^)®[^W^(-A,7,0m)«i*(7r)]7 



xl a 0^ fc/i (7r)i? jV (7r)^(0 m ) 

xl a V d U z (-^ m )R kfl (7v)R ju (7v). 
with k : j = l or 2, (fc ^ j), and /i, v = x or 2/, (p ^ z/). 



3. Phase measurement: Quadrature detection in 
NMR 



After NMR interferometry, the state of the quantum 
register is ^(|0) a +e i(/ ^+Ai) |l) )|*s)i2- A relative phase 
shift Pg + Pd is created between the states |0) a and |l) a of 
the auxiliary qubit. It can be obtained when we measure 
the NMR signal of the auxiliary qubit: 

(O = - i° a y ) = \[™s{P g + fa) + ism(fa + fa)]. 

The quadrature detection in NMR serves as a phase 
sensitive demodulation technique, i.e., the complex de- 
modulated signal is separated into two components (the 
real part RE ex cos(/3 g + (3d) and the imaginary part 



c 

.Sf 

as 

z 



10} + |i> 



3 H> 



2 H> 





NMR frequency 



FIG. 9: Phase measurement by quadrature detection in NMR. 



IM oc sm(/3 g + Pd)) which are 90 out of phase with each 
other. Thus, the phase angle j3 g + /3d of the signal can 
be determined by tan -1 (1 M / RE) . Quadrature detection 
combined with Fourier analysis thus gives all the neces- 
sary information on the magnetic resonance signal com- 
ponents i.e., amplitude, phase and frequency [28]. Taking 
the input state of 775 (|0) a + |l) a )|^#)i2 as the reference 
spectrum, we measured the relative phase information by 
the phase of the Fourier-transformed spectrum. Fig. [9] 
shows a simple example for the phase measurement by 
the Fourier-transformed spectra. 



Experimental results and data analysis 



Fig. 10 shows the experimental NMR spectra for a 
set of experiments H(X : j) with varying magnetic field 
strength A (0 — >• 1.7) and a constant anisotropy parame- 
ter 7 = 0.5. The parameter A was varied by a hyperbolic 
sine function [30 , but avoiding the level crossing points. 
)T From these spectra, we measured the phase shifts accu- 
mulated by the two trajectories C and C listed in Table 
I. The pure GP was then obtained by summing the two 
phase shifts from C and C: f3 g = (f3 t (C) + /3 t (C))/2. A 
set of the experimental spectra are shown in Fig. 3 in 
the paper. 

Our experimental errors of the geometric phases are 
less than 3°. These errors result from the imperfection 
of the initial ground state by ASP, the diabatic effect, 
and other experimental imperfections such as the inho- 
mogeneity of the radio frequency field and the static mag- 
netic field, and the imperfect calibration of the radio fre- 
quency pulses. The decoherence from spin relaxation was 
small, since the total experimental time of less than 90 
ms was short compared to the shortest relaxation time 
of 1.0 s. 

The error contributed by the imperfection of the ini- 
tial ground state of ASP can be evaluated by the use of 
the measured input density matrices p^f after ASP, e.g, 
shown in Fig. [7] Therefore we have the input state 



|0) a (0|<8>p: 



exp 



10 




FIG. 10: Experimental 13 C spectra for the XY Hamiltonian 
7-L(\, 7) with 7 = 0.5 and the variable A from to 1.7. The red 
ones are the spectra of the initial states prepared by ASP as 
the references. The upper plot is the spectra for the trajectory 
C when A scans while the lower plot is for C. We extract 
the phase information from these spectra using quadrature 
detection. 



TABLE I: The extracted phase values from the NMR signals 
for 7 = 0.5. f3^ p and f3^ p denote the phases obtained by the 
NMR interferometry with the related adiabatic evolution tra- 

(3 exp +(3 e - xp 

jectories C and C. f3g Xp = — — 2 c is the geometric phase. 
The corresponding experiment spectra are showed above. 



ExpNo. 


A 


Pc P (°) 


p% p (°) 




P* h (° 








170.6 


-173 


-1.2 





1 


0.327 


174 


-171.4 


1.3 





2 


0.5312 


173.4 


-172 


0.7 





3 


0.6592 


173.8 


-174.2 


-0.3 





4 


0.74 


168 


-170 


-1 





5 


0.7922 


172.2 


-171 


0.6 





6 


0.8275 


174.6 


-176.2 


-0.8 





7 


0.8541 


170.4 


-171.6 


-0.6 





8 


0.878 


-62.6 


114.2 


25.8 


23.6 


9 


0.9045 


-60 


107.6 


23.8 


22.5 


10 


0.9399 


-60.2 


106.4 


23.1 


21.1 


11 


0.992 


-69.8 


111.4 


20.8 


19.3 


12 


1.0729 


-70.4 


108.6 


19.1 


16.8 


13 


1.2009 


-79.6 


109 


14.7 


13.8 


14 


1.4051 


-75.6 


95.4 


9.9 


10.4 


15 


1.7321 


-81.1 


96.6 


7.7 


7.1 



Then we input the state to an ideal NMR interferom- 
eter, i.e, a perfectly implemented Hadamard gate and 
controlled evolution process simulated on a classical com- 
puter to get the theoretical output: 

if pTJ &Ut) \ 

Pf 2 ^ U a (T) P rj U a {T)p?JUt{T) ) ' 

where a (a = C or C). As a result, the measurement on 
the auxiliary qubit by the quadrature detection gives 

(p f a-) = l -Tr[U a {T)pZ^ 

Thus we achieved the simulated phases f3 s ^ m = 
arg((p f a-)) = arg[Tr(U a (T)p™f)}. Thus the simulated 
geometric phase starting from the experimental initial 

osim I osim 

state is (3 s g irn = — — 2 c . We found that the errors 
contributed by the imperfection of the prepared initial 
state is about 1°. 



References 

[24] H. F. Trotter, Pacic J. Math. 8, 887 (1958). 

[25] I. L. Chuang et al., Proc. R. Soc. A 454, 447 
(1998); D. G. Cory, A. F. Fahmy, and T. F. Havel, 
Proc. Natl. Acad. Sci. U.S.A. 94, 1634 (1997). 

[26] G. M. Leskowitz and L. J. Mueller, Phys. Rev. A 
69, 052302 (2004). 

[27] X. Peng and D. Suter, Front. Phys. China 5, 1-25 
(2010). 

[28] R. R. Ernst, G. Bodenhausen, and A. Wokaun. 
Principles of Nuclear Magnetic Resonance in One 
and Two Dimensions. Oxford University Press, 
Oxford, 1994. 

[29] M. V. Berry, Proc. R. Soc. Lond. A 392, 45 
(1984). 

[30] X. Peng et al, Phys. Rev. A 71, 012307 (2005). 



